Charge transport-mediated recruitment of DNA repair enzymes 
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Damaged or mismatched bases in DNA can be repaired by Base Excision Repair (BER) enzymes 
that replace the defective base. Although the detailed molecular structures of many BER enzymes 
are known, how they colocalize to lesions remains unclear. One hypothesis involves charge transport 
(CT) along DNA [Yavin, et ai, PNAS, 102, 3546, (2005)]. In this CT mechanism, electrons are 
released by recently adsorbed BER enzymes and travel along the DNA. The electrons can scatter 
(by heterogeneities along the DNA) back to the enzyme, destabilizing and knocking it off the DNA, 
or, they can be absorbed by nearby lesions and guanine radicals. We develop a stochastic model to 
describe the electron dynamics, and compute probabilities of electron capture by guanine radicals 
and repair enzymes. We also calculate first passage times of electron return, and ensemble-average 
these results over guanine radical distributions. Our statistical results provide the rules that enable 
us to perform implicit-electron Monte-Carlo simulations of repair enzyme binding and redistribution 
near lesions. When lesions are electron absorbing, we show that the CT mechanism suppresses 
wasteful buildup of enzymes along intact portions of the DNA, maximizing enzyme concentration 
near lesions. 

PACS numbers: 87.15.H,82.39.Pj,05.10.Gg,05.40.-a 



I. INTRODUCTION 

The genomes of all living organisms arc constantly un- 
der attack by mutagenic agents such as reactive oxygen 
species and ionizing radiation. Such processes can dam- 
age bases giving rise to localized lesions in the DNA^i^ 
that can lead to harmful mutations and diseases such as 
cancer. For example, guanine residues can be oxidized, 
generating a radical called 7,8-dihydro-8-oxoguanine, or 
oxoG for short. Unlike the non-oxidized form, this radi- 
cal can pair with both cytosine and adenine, uhimately 
giving rise to GC — > TA transversion mutationsi upon 
muhiple replications. Lesions can also arise through alky- 
lation, hydration and deaminatioufi 

One defense mechanism against these mutation pro- 
cesses is the Base Excision Repair (BER) pathway. BER 
enzymes recognize and undo damage to DNA by adsorb- 
ing onto the sugar-phosphate backbone, locating the le- 
sion and excising it. The biomechanical functions of re- 
pair enzymes have been well established and their 3D 
structures are known in great detaili^ There are four 
main types of BER enzyme: DNA glycosylases, AP- 
endo nucleases, DNA polymerases and DNA ligases. Each 
of these enzymes has a different role in the BER family. 
For example, DNA glycosylases initiate the repair path- 
way, detecting and recognizing distinct forms of DNA 
damage while the endonucleases are responsible for cleav- 
ing the sugar-phosphate backbone. Together, these en- 
zymes maintain the overall integrity of DNA, generally 
ensuring that miscoded proteins are kept to a minimum. 

The problem of how a BER enzyme locates a lesion 
on DNA is a specific example of how enzymes find lo- 
calized targets. The DNA of E. coli contains about 10^ 



base pairs. If we assume that BER enzymes find lesions 
through a pure ID diffusive "sliding" process with diffu- 
sion constant D base pairs^/s, the search time is roughly 
10^^/D. Estimating D to be 5 x 10^ base pairs^/s, the 
value for a human DNA glycosylase;^ we obtain a search 
time of about 2 x lO^s w 2 days, much longer than even 
the reproductive period of E. coli. Therefore, it is likely 
that other mechanisms are responsible for DNA target 
location. 

In 1970, Riggs et alXi^ measured the association 
rate of the Lad repressor protein to its target on 

This was puzzling 



DNA to be about lO^^M'^s 



because the theoretical upper limit for the association 
rate of a Lad enzyme diffusing in 3D is predicted (via 
the Debye-Smoluchowski formula) to be about 2 or- 
ders of magnitude less. This fundamental biophysi- 
cal problem was studied in the seminal work of von 
Hippel and co-workers2ii£iiiii2 and the "faster-than- 
diffusion" search of targets on DNA has received re- 
cent attention] -'^'^i-'^^'^^i-'^^i^^i^^ Facilitated diffusion is one 
mechanisroiiii^ii^ii^ proposed to explain the accelerated 
search. Instead of diffusing directly to their target, the 
searching enzymes can spend part of their time attached 
to the DNA and perform a ID random walk along part of 
the strand. If the enzyme is able to spend 50% of its time 
on the DNA and 50% of its time diffusing in 3D, and the 
diffusion constants in ID and 3D arc comparable, the 
association rate is predicted to increase by as much as 
lOOf^ bringing it in line with the experiments in Riggs 
et al^^ However, other authors have shown that (i) typ- 
ical enzymes are highly associated with DNA, spending 
over 99.999% of their time on the strand^ and (ii) the 
diffusion constant in ID can be 1000 times smaller than 
in 3D)^ resulting in a negligible reduction of the search 
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FIG. 1: Redox mechanism for repair enzyme interaction based 
on the papers by Yavin, et alA and Boon, et al^ (a) A MutY 
in the 2+ state (solid hexagon) adsorbs and oxidizes to the 
3+ state (empty hexagon) by releasing an electron along the 
DNA. The electron is emitted to the left or right of the en- 
zyme with equal probability. Guanine radicals (circumscribed 
dots) can absorb electrons and prevent oxidation of nearby 
adsorbed enzymes, (b) A MutY repair enzyme in the 3+ 
state absorbs an electron and is reduced, causing it to desorb. 
(c) Lesions also prevent passage of electrons, either through 
electron absorption or reflection. In our analysis, lesions act 
differently from oxoG radicals in that they can continuously 
absorb electrons. 



time. Hence, facilitated diffusion in its basic form is not 
adequate to explain the fast reaction rates observed. Ex- 
tensions to the facilitated diffusion theory can incorpo- 
rate finite enzyme concentrations "antenna" effects re- 
sulting from the conformation of the DNA^ii fast inter- 
segment transfers of the protein^i^ specifrc/non-specific 
protein-DNA interactions and directed DNA sliding 3i 

Although BER enzymes may colocalize to lesions 
by exploiting the facilitated diffusion mechanisms cited 
above, other mechanisms arc likely required for efficient 
and timely recruitment to lesions. A charge-transport 
(CT) mechanism has been recently proposed as a possi- 
ble basis for efficient scanning by MutY, a type of DNA 
glycosylase4i^ MutY is known to contain an iron-sulfur 
cluster which plays a key role in the CT mechanism. The 
cluster can take one of two forms: [4Fe-4S]^+ and [4Fe- 
4S]'^+. When MutY is in solution, the cluster is in the 2+ 
state and is resistant to oxidation. However, upon bind- 
ing to DNA, the cluster potential is shifted, making the 
3-1- state more accessible. The result is that after bind- 
ing, MutY-[4Fe-4S]^"'" is easily oxidized and releases an 
electron along the DNA, as shown in Fig[T{a). It should 
be noted that the 3-|- state of MutY has a binding affinity 
that is about 4 orders of magnitude larger than that of 
the 2+ stated Therefore MutY-[4Fe-4S]^+ spends most 
of its time in solution whereas MutY-[4Fe-4S]'^"'' exists 



primarily adsorbed onto DNA. 

Although controversial about 15 years ago, long 
range electron transport in DNA is now a well accepted 
phenomenonj ^^i^^ Experiments indicate that charge 
transport can occur over 40A (about 12 base pairs) 
in less than a nanosecon d^^i^^ and the influence of 
DNA strand crossovers on CT is generally smalli^ 
Although electron dynamics along DNA is in general 
very complicated, some aspects of the process are now 
understood. For example, both guanine and adenine 
can act as carriers of positive charge; in analogy with 
semiconductors, oxidized DNA can transport charge via 
the transfer of holes from base to base. 

Quantifying how BER enzymes adsorb to DNA and 
how they are recruited to lesions has so far been restricted 
to simple scaling argumentsi^ In this paper, in order to 
explore the implications of DNA target selection solely 
by CT, we assume that adsorbed MutY BER enzymes 
do not slide along the DNA. However, upon first attach- 
ment to DNA, the enzyme will emit an electron that 
propagates along the strand in a random direction and 
its cluster wiU go from the [4Fe-4S]2+ to the [4Fe-4S]^+ 
state. Should this electron become absorbed by another 
MutY-[4Fe-4S]^+ enzyme further along the DNA, the 3+ 
form is reduced and desorbs (Fig[TJb)). If the electron 
back-scatters and returns to the original MutY, it self- 
desorbs. Although the model proposed in this paper is 
intended to specifically describe the colocalization and 
redistribution of MutY through the redox reaction of its 
iron-sulfur cluster, many BER enzymes, in fact, contain 
such a cluster, e.g. endonuclease HI. Therefore, we think 
that our model may be more general and could also de- 
scribe the binding kinetics of other enzymes. 

Since unbiased stochastic motion in ID always leads to 
return of the electron;^ in the absence of any other elec- 
tron absorbers on the DNA, a MutY BER enzyme that 
is deposited will eventually self-desorb with probability 
1. However, BER enzymes can be recruited to DNA by 
preexisting electron absorbers. These are typically gua- 
nine radicals ("oxoG") and other lesions, indicated in 
Fig. [TJc) by circumscribed dots and filled diamonds, re- 
spectively. It has been suggested that oxoG plays an 
important role in the seeding of MutY onto DNAj^ The 
oxoG radicals, like adsorbed enzymes, are able to absorb 
electrons, preventing them from returning and desorbing 
BER enzymes that originally released them. Therefore, 
the oxoG radical in Fig. [TJa) can absorb one left-moving 
electron and prevent it from back-scattering and desorb- 
ing the right-most enzyme. Upon reduction, oxoG rad- 
icals convert to irormal guanine bases, no longer absorb 
electrons, aird no longer take part in the CT mechanism. 

Other lesions do not simply annihilate by absorbing 
electrons; rather, they require the physical presence of 
BER enzymes to excise them. These lesions may recruit 
smaller, more abundant proteins from solution that per- 
mit multiple electron absorption. Another possibility is 
that the lesions reflect electrons. Both cases are shown 
in Fig. [ijc). Therefore, our basic model consists of right 
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and left-moving electrons, guanine radicals, oxidized and 
reduced forms of BER enzymes, and lesions on the DNA 
strand. Newly adsorbed BER enzymes instantly release 
electrons (right or left-moving), while oxoG radicals, le- 
sions, and oxidized BER enzymes absorb electrons and 
prevent their passage. 

In this paper, we model the adsorption, desorption and 
redistribution of repair enzymes using the redox mecha- 
nism shown in Fig. [T] Wc first derive some exact results 
in the absence of any lesions; in particular, enzyme ad- 
sorption probabilities and the time taken for returning 
electrons to induce enzyme desorption. These results en- 
able us to define rules for Monte-Carlo simulations of the 
dynamics of multiple enzymes. For electron absorbing 
lesions, simulations show that if enzymes are deposited 
onto a DNA at a rate that is slow compared to the elec- 
tron dynamics, the distance between a lesion and the 
closest enzyme scales as 0(n~^/^) for large n, while total 
number of enzymes adsorbed between two lesions scales 
as 0(77.^/"^). However, because of the CT mechanism, 
this accumulation is not uniform along the DNA and 
the maximum enzyme density always occurs at lesions. 
Hence for electron-absorbing lesions, the CT mechanism 
concentrates enzymes to damaged bases in DNA, consis- 
tent with the qualitative predictions in Yavin, et alA and 
Boon, et al^ 

The outline of this paper is as follows. In the next 
section, we develop a model for the electron dynamics 
based on the stochastic Broadwell model;^'^'^'^ Pairs 
of guanine radicals, BER enzymes or lesions define the 
boundary of a segment (a "gap" ) over which an electron 
can propagate. Section HI contains our results. In Sec- 
tion llll Al we we derive enzyme sticking probabilities and 
the time taken for returning electrons to desorb the en- 
zymes that originally emitted them. In particular, we 
derive the MutY desorption rate in terms of the electron 
scattering (flip rate) and the electron speed. In Section 
IIIIBl we perform implicit-electron Monte-Carlo simula- 
tions to study the redistribution and accumulation of en- 
zymes between two fixed lesions on the DNA. Finally, in 
Section IV, wc discuss facilitated recruitment of enzymes 
to lesions in the context of the CT hypothesis, as well as 
the biological advantages and disadvantages of the pro- 
posed CT mechanism. 



II. STOCHASTIC CHARGE TRANSPORT 
MODEL 

A. One-sided Broadwell problem 

In analogy with Bicout's analysis for the unrelated 
problem of microtubule growth dynamics,— we now 
present similar equations for the dynamics of electrons as- 
sociated with repair enzymes. Consider Fig. [2l^a): oxoG 
guanine radicals with density p are distributed randomly 
along an infinite strand of DNA. A single repair enzyme 
initially attaches to the DNA at a random position, in 



Symbol 


Definition 


units 


p+ 


Probability density of rightward 


1/L 




cicC LI (Jli 


P- 


J. i.\J KjOL/t-jLLi. Vi y V-i.t'llOl Ij V w± 1C± L VV Cl± \JL 

electron 


1/L 


X 


Position along DNA 


L 


Xq 


Position of electron release 


L 


T 


Time 


T 


P 


Density of oxoG guanine radicals on DNA 


1/L 


L 


Distance between two oxoGs/enzymes 


L 


F 


electron flip rate 


1/T 


V 


Electron speed 


L/T 


AI 


Electron decay rate 


1/T 


h 


Deposition rate of enzymes 


1/(L . T) 



TABLE I; Table of dimensional variables and parameters. 
The analysis performed assumes A/ = 0. L represents length 
and T represents time. 
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Math defn. 


Descriptive definition 


Q± 


P±/P 


Rightward/ Leftward electron 
probability density 


X 


pX 


Coordinate along DNA 




pXo 


Position of electron release 


t 


pVT 


Time 




pL 


"Gap size" : distance between two 


£ 


oxoGs / enzymes /lesions 


f 


F/pV 


Electron fiip rate 




M/pV 


Electron decay rate 






Position of enzyme adsorption 






Enzyme- lesion /enzyme-enzyme 




distance (see Fig. [SI 



TABLE II; Definitions of dimensionless symbols in terms of 
the dimensional quantities in Table HI 



between two electron absorbing oxoGs. The enzyme im- 
mediately emits an electron along the DNA to the left or 
right with equal probability. The electron can only move 
with speed V , in the positive or negative X— directions, 
executing random flips between the two directions with 
rates F. Furthermore, emitted electrons can be annihi- 
lated with rate M through nonspecific interactions with 
random electron absorbers diffusing in the bulk. 

In general, two steps arc required for a MutY enzyme 
to bind to DNA. First, when MutY-[4Fe-4S]^+ is in con- 
tact with the DNA, it has to undergo oxidation by re- 
leasing an electron. The oxidized form of the enzyme 
binds more strongly to DNA. Second, the released elec- 
tron must be absorbed by some particle other than the 
enzyme (an oxoG, an already adsorbed MutY or a lesion) 
to prevent it from returning and reducing the enzyme. 
This allows the enzyme binding to become "permanent" . 
Therefore the net binding probability depends on (i) the 
probability of electron r eleaseby MutY-[4Fe-4S]2+ (when 



4 



(a) 


^ e~ 




















-zy2 



(c) 



FIG. 2: (a) A repair enzyme (hexagon) adsorbs onto a DNA 
which is initially populated by guanine radicals (circled dots) 
with a density p. An electron is emitted to the left or right 
with equal probability. The emitted electron has flip rate F, 
speed V , and decay rate M. (b) The one-sided Broadwell 
problem. An electron is emitted from X = with probabil- 
ity 1 toward a guanine radical sX, X = L. (c) The two-sided 
Broadwell problem. An enzyme is deposited between two gua- 
nine radicals which are a distance L apart. Immediately after 
landing inside this segment, an electron is emitted to the left 
or right with equal probability. 



in contact with the DNA) and (ii) how far neighboring 
electron absorbers are from the adsorbed MutY. In this 
paper, we assume that when enzymes adsorb onto the 
DNA, they always oxidize, releasing an electron with 
probability 1; in Fig. HJb), the electron is released to 
the right with probabihty 1 and in Fig. [D^c), the elec- 
tron is released to the left or right with probability 1/2. 
In principle, an enzyme can attach to and then immedi- 
ately detach from the DNA without releasing its electron, 
but assuming the electron release rate is large, we neglect 
this process. The adsorption probabilities we derive later 
in this section will depend only on the gap size L and the 
parameters for electron motion. 

Finally, we assume that the DNA is immersed in an in- 
finite reservoir of enzymes which is kept at a fixed chemi- 
cal potential. The rate of deposition of enzymes onto the 
DNA is assumed to be constant. A deposited enzyme can 
either adsorb by having its released electron captured by 
neighboring electron absorbers or it can desorb due to its 
electron returning. 

To build our full solution, wc first derive exact analyti- 
cal expressions for the "one-sided" problem shown in Fig. 
[2jb) which consists of an enzyme at X = and a guanine 
radical at X ~ L. At time T = 0, an electron is emitted 
from a position Xq > (subsequently, wc will take the 
limit Xq 0) with speed V in the positive X-direction. 
For the one-sided problem, the electron is emitted only 
to the right. The probability that the electron is at a po- 
sition between X and X + dX at time T, and moving to 
the right with velocity V is denoted P+ {X, T) . Similarly, 
P_ (X, T) denotes the probability density of an electron 
moving with speed V in the negative X-direction. The 
electron can flip directions by scattering from inhomo- 
geneities and thermally excited conformational variations 



along the DNAj^i^i We model this flipping process as 
a spatially homogeneous process occurring with constant 
rate F, independent of any structure along the DNA such 
as base pair sequence. 

The evolution equations for the probability densities 
P±(X,T) are 



dP+ 
dP- 



-V 



dP. 



+ 



dX 

dP- 
'dX ' 



FP+ + FP- 



MP+ 



(1) 



FP+ - FP- - MP- 



where < X < L. Eqs. ^ describe the probability 
density of electrons being advected to the right and left. 
The flipping of the electron motion is represented through 
F and couples the equations for P+ and P-. Further- 
more, the densities decay in time with an annihilation 
rate M. Electrons can be annihilated by being absorbed 
by other proteins (besides BER enzymes) in solution. If 
these proteins adsorb onto the DNA, absorb an electron 
and desorb back into solution, an electron is permanently 
removed from the DNA. 

The boundary conditions and initial conditions are 



p+(o,r) = p_(L,r) = 0, 
p+ix,o) = S{X 
P-{X,0) = 0. 



(2) 
(3) 
(4) 



The boundary conditions ([2]) arise because the enzyme at 
X = and the oxoG at X = L (see Fig. [DJb)) are both 
perfect electron absorbers. When Xq 0, the initial 
condition ([3]) reflects the fact that an electron is released 
to the right from the enzyme at A" = 0. Initially, there 
arc no leftward traveling electrons in Fig. [IJb), justifying 
Eq. ([4]). All variable and parameters are listed in Tables 

m 

We now define dimensionless independent variables 
through the guanine radical density p and the rightward 
electron travel time l/{pV): 



X = pX, t = pVT, 
so that Eqs. ([T]) can be written in the form 

Q+{x,t) 

Q-{x,t) 
where Q± = P±/p and 



(5) 



(6) 



f 



(7) 



and < X < £ = pL. In Eq. 

F_ 

^ ^ pV ^~ pV 



(8) 
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is the dimensionlcss flipping rate and electron decay rate. 
The boundary and initial conditions (O, ©J (HI) become 



-(o,i) = g_(£,t) =0, 

Q+{x,0) ^S{x 
Q-{x,0) =0, 



xo) 



(9) 



where xq = pXq. In the physical problem, an electron is 
released from the enzyme as soon as it initially attaches 
to the DNA. Therefore, we solve Eqs. ([6]) with ^ and 
dn]) taking the limit xq ^ (for details, see Appendix [A)) . 
The dimensionless variables are tabulated and defined in 
Table HIl Henceforth all of our results and analyses will 
be presented for /i = 0. 

The probability of the enzyme in Fig. ^h) sclf- 
desorbing before time t is given by Q_(0,t')dt' where 
Q- can be found by taking the inverse Laplace Transform 
of Eq. (jA4p in Appendix \X\ Therefore, the enzyme des- 
orption and sticking probabilities for the one-sided prob- 
lem are 



1 + fi 
respectively. 



and 1 — 



1 



l + /f l + fV 



(10) 



B. Two-sided Broadwell problem 

Now consider the two-sided problem depicted in Fig. 
[5]Jc). A repair enzyme lands at position ^ between two 
oxoG guanine radicals that are a distance I apart. The 
solution to the full problem can be found by splitting it 
into two subproblems and using our results from Section 
III Al Instead of solving for the densities on [0,£], we can 
solve for Q± separately on [S_,£/2] (with the enzyme ini- 
tially deposited at ^ and the guanine radical at ^/2), on 

— ^/2] (with the enzyme at ^ and the guanine radical 
at — ^/2) and combine the results. The enzyme desorp- 
tion and adsorption probabilities pO)) extend straightfor- 
wardly: 



dcsorb 



1 

2 

nadsorb(^, ^) = 9 



1 



l + /(f-c) l + /(l+o 

/(II) 

A plot of the sticking probability Iladsorb for different 
values of / and for two different gap sizes is shown in 
Fig. [3l For a fixed gap size, and sufficiently large / (cor- 
responding to a diffusive electron motion), permanent 
BER enzyme adsorption is less likely to occur near the 
center of the gap because absorption of the electron by 
guanine radicals is less likely to occur. The permanent 
adsorption or sticking probability is more uniform when 
/ is small (corresponding to a ballistic electron motion): 
whether the oxoG radical is close or far away from the 



Symbol 


Descriptive definition 


See Eq. 


nadsorb(C5 


Enzyme adsorption probability 


m 


-^-^QGSOrD \S 5 / 


Enzyme desorption probability 


(fTTI) 


Iladsorb (i*) 


Enzyme adsorption prob. averaged 
over landing position ^ 


m 


ndesorb(-^) 


Enzyme desorption prob. averaged 
over landing position ^ 


(Uni) 




Enzyme adsorption prob. averaged 
over landing posn. ^ and gap size i 




(Iladsorb ) 


(|14|16|) 


u 


Random variable for conditional 
return time of electron 


m 




Mean conditional return time 


m 


(MCRT) of an electron 




MCRT of an electron averaged over 
landing position ^ 


m 


(rr) 


MCRT of an electron averaged over 
landing posn. ^ and gap size £ 


m 



TABLE III: Derived adsorption/desorption probabilities, 
electron return times and related quantities. 



enzyme makes little difference to the adsorption probabil- 
ity. Finally, for fixed /, increasing the gap size decreases 
the adsorption probability because guanine annihilation 
by the electron is less likely to occur. The diffusive and 
ballistic behaviors of the Broadwell model are derived in 
Appendix |B] 




FIG. 3: Dependence of enzyme sticking probability, Iladsorb 
(see Eq. (HU), on dimensionless flip rate / and landing po- 
sition —£/2 < ^ < £/2, for the deposition of a single enzyme 
into a gap (the segment of DNA between two guanine radi- 
cals) of size £. (a) ^ = 1 with radicals located at ±0.5 (b) 
£ = 2 with radicals located at ±1. 



III. RESULTS AND DISCUSSION 

A. Statistics of repair enzymes away from lesions 

In this section, we present and discuss deposition 
statistics that are valid far away from lesions. First, 
using Eq. (jlip . we average over the landing position 
^ to calculate mean sticking/adsorption probabilities of 
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repair enzymes that are deposited between two guanine 
radicals that are a distance I apart. The inter-radical 
distances ("gaps") in DNA will, in general, be randomly 
distributed. Therefore we ensemble-average our results 
over the distribution that i is expected to obey. Sec- 
ond, we find the mean return times of electrons, i.e. the 
time taken for a deposited enzyme to be desorbcd by its 
own electron, providing it dcsorbs. Again, our results arc 
ensemble-averaged over randomly distributed gap sizes. 
The quantities we shall compute and analyze in this sec- 
tion are listed in Table Hill 

All the results presented are for adiabatic depositions. 
A deposition is adiabatic if the inter-deposition time is 
much larger than the time scale of the electron dynamics. 
In other words, for every enzyme deposited, its released 
electron completes its motion before the deposition of the 
next enzyme. At any given time, there is at most one 
traveling electron on the DNA. For details, sec Appendix 



1. Repair enzyme sticking probability 

One quantity of interest is the probability that any 
given repair enzyme that lands on the DNA will not be 
kicked off by its own electron, and will remain adsorbed. 
Enzyme sticking relics on efficient capture of the released 
electron by neighboring electron absorbers (guanine rad- 
icals and adsorbed enzymes). Intuitively, one would ex- 
pect that a greater density of absorbers with smaller gaps 
would result in a more efficient capture of enzymes. 

For a single repair enzyme deposited onto the DNA, 
landing at a position —i/2 < C < (see Fig. [H^c)) 
inside a gap of length £, centered about a; = 0, the prob- 
ability of it remaining on the DNA is given by Iladsorb in 
Eq. (fTT|) . This quantity can be averaged over all possible 
deposition positions ^ within the gap to obtain 



adS' 



1 

2e 



i/2 



1 



1 



■ tanh 



/(I 



i + /(f + ?) 



V2 + /^. 

(12) 

This result is plotted in Fig. g] (dashed hue). Eq. 
gives the sticking probability of a repair enzyme newly 
deposited between two electron absorbers separated by 
£, uniformly averaged over its deposition position within 
the gap. 

We now average over the gap length distribution to 
compute the sticking probability for deposited enzymes 
that land anywhere along the entire DNA strand. For 
an infinite, lesion free DNA, depositing an enzyme will, 
in general, change the local guanine and enzyme distri- 
bution. Hence, the sticking probabilities will also change 
with each successive deposition, making the calculation 
difficult in the context of the Broadwell model. However 
it is possible to calculate the sticking probability for a 



given gap distribution. In special cases where this distri- 
bution is known or simple to calculate, we can compute 
the efficiency of enzyme recruitment onto the DNA. 

Consider the case of a DNA with a discrete distribu- 
tion of gaps £i,£2,£3, ■■■■ Suppose that on the DNA, a 
fraction of the gaps have size £j. Now consider many 
realizations of a single enzyme deposited onto this DNA. 
The fraction of enzymes that lands in gaps of size £j is 
/ ^j=i 't'j^j ^^^d the fraction of these that stays ad- 
sorbed, using Eq. is 



2 ' 



; tanh" 



-1 (Jl^\ 



(13) 



The fraction of enzymes that stays adsorbed (in any gap) 
is obtained by summing over j. In the continuum limit, 
£j £, (pj (l){£)d£, where £ is the continuous gap length 
and (j){£) is the probability distribution function (PDF) 
for £. We obtain 



f{e) 



(j){£) tanh" 



2 + f£ 



d£. (14) 



Note that (Hadsorb) 7^ /o°° </'(^)nadsorb(^), the resuh that 
one might expect by naively averaging Eq. (|12p over the 
gap distribution. 

Since (j){£) depends on the number of enzymes de- 
posited, it is time dependent. In principle, one could 
calculate how (j)(£) changes as enzymes are adiabatically 
deposited. The corresponding evolution of the sticking 
probability is then given by Eq. (|14p . One possible way 
of finding how (j){£) evolves is to use a mean field the- 
ory for the particle distributions, but we leave this as the 
subject of a future investigation. 

In the special case where one enzyme is deposited onto 
a DNA that only has guanine radicals, we can calcu- 
late (j){£) and hence (Hadsorb) explicitly. If the guanine 
radicals have a number density p, then the gap lengths, 
on average, are 1/p, which corresponds to a imit dimen- 
sionless gap size (see Eq. (O). Hence {£) ~ 1 and the 
dimcnsionless gap sizes, Y, are exponentially distributed 
(see Appendix [P]) according to 



Prob(£ < gap size < £ + d£) = e ^d£, 



(15) 



so we set (p{£) — e 
([13]), we obtain 



Substituting this result into Eq. 



(n, 



adsorb/ 



e^/^Ei(l//) 
/ 



(16) 



where Ei(a;) = ^l~dt is the exponential integral. This 
analytic result is plotted in Fig. |4] (solid line) and is con- 
firmed by Monte-Carlo simulations (circles). The sticking 
probability increases when either the electron-absorber 
density p increases, the electron velocity V increases or 
the flip rate F decreases. 

Equation ()16p is valid only when the number of en- 
zymes that have stuck is much less than the initial num- 
ber of oxoG radicals. In this limit, the distribution of 



7 



5 -"^ 10 15 
1 1 I 




I I 

5 10 15 

flipping rate f=F/(p V) 



enzyme sticking probability inside an inter-radical gap 
of unit length, whereas (Iladsorb) (Eq. ([TB)) ) gives the 
enzyme sticking probabihty averaged over exponentially 
distributed inter-radical gaps lengths, but with unit 
mean. Intuitively, one would expect the boundaries 
defining the smaller gaps to be more efficient at seques- 
tering electrons than those associated with larger gaps. 
However, the enhanced electron trapping by smaller gaps, 
leading to otherwise increased sticking probabilities is 
compensated by a higher deposition flux into larger gaps 
(large gaps collect more enzymes than small gaps). The 
net result of averaging over exponentially distributed gap 
sizes is for the larger gaps to dominate and lower the 
overall gap-averaged sticking probability. This is shown 
in Fig. m where for aU values of /, (Iladsorb) < Iladsorb 
when i = 1. 



FIG. 4: Enzyme sticking probabilities as a function of di- 
mensionless flip rate /. The simulation data was obtained by 
performing single depositions onto a DNA of length 100 {i.e., 
with physical length lOQ/p). The fraction of enzymes that 
remain on the DNA after performing 10''' trials was recorded. 
Increasing the DNA length did not significantly affect the 
simulation results. 

gap lengths will remain approximately exponential. For 
the human genome of ^ 10^ base pairs, there are approx- 
imately 10'' oxoGs present at any given timei^ In this 
case, we expect that Eq. should be fairly accurate 
for about the first dozen depositions. 

Note that fladsorb (Eq. (dH)) with £ ^ 1 gives the 



2. Mean conditional return time of electrons 

We now find the mean time that a BER enzyme stays 
on the DNA after its initial deposition, conditioned on 
its own electron returning and knocking the enzyme off. 
This quantity allows us to estimate a rate of desorption 
that can be used in more coarse-grained, higher level de- 
scriptions of the CT mechanism. 

Consider depositing an enzyme into a gap of size £ at 
a position ^ satisfying —£/2 < ^ < £/2. The probability 
that the electron ("e~") returns in a time < t, given 
that it returns is. 



w I - ^ Prob(tr < t) 

Prob(tr < t I e returns) — 



Prob(e returns) 
Proh{tr < t \ e~ shoots right) + Prob(tr < t \ shoots left) 



Prob(e returns | e shoots right) -I- Prob(e returns | e shoots left) 

(17) 



i /o Q_ (x = 0, t'; 0, £/2 - ^dt' + \ /J Q_ (x = 0, i'; 0, + i)dt' 



i Q_ (x = 0, t'; 0, £12 - e)dt' -f- \ Q_ {x = 0, t'; 0, £/2 + i)dt' ' 



In Eq. pT|) , Q- (a;, t] xq , £) is the leftward electron density 
at position < a; < ^ at time t given that the electron 
was released from a; = at i = (see Fig. [Ub) for the 
xq = case). This density comes from solving Eqs. © 
and ([7]) along with the conditions 

The mean conditional electron return time can then 
be computed from 

/■°° d 

Tr{i;£J) = t-^ Vroh{tr < t\e^ returns) dt. (18) 
Jq 

Using Eq. dni), Tr{^;£,f) in Eq. ^ can be found in 
terms of the Laplace-transformed density Q±{x, s) which 



is given in Eq. (jA4p of Appendix Upon averaging 
''r(C; £, /) over the initial landing positions ^, we obtain 



f^(^, /) = A /^ + ^^ tanh-i 



3/ ^f£{2 + fi) 



fi 



+ ^[/^-l--^log(l + /^)]. 



(19) 



We plot fr{£, /), and validate Eq. p9|) using MC simula- 
tions in Fig. EJa). Finally, we further ensemble- average 
fr over gap lengths £. Consider many realizations of the 
deposition of a single enzyme onto an infinite DNA with 
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FIG. 5: (a) For adiabatic enzyme depositions into a gap of 
size £, fr, the mean conditional return time of an electron, 
averaged over the enzyme landing position ^, is recorded for 
different gap sizes and dimensionless flip rates /. The symbols 
represent data from Monte Carlo simulations and the solid 
line represents the analytic expression from Eq. (|19|l . (b) For 
random, uniform, adiabatic enzyme depositions onto a DNA 
with randomly and uniformly distributed guanine radicals, 
(fr) as predicted by Eq. is plotted as a function of the 

dimensionless flip rate /. 

oxoGs whose gaps are exponentially distributed. The av- 
erage time that the enzyme stays adsorbed, given that its 
electron eventually returns to knock it off, is (fv). The 
calculation of (fv) is similar to that of (Iladsorb) described 
in Section IIII A 11 but modified to account for the fact 
that the number of enzymes that self-dcsorb (i.e., the 
number of return times that are finite) depends on £. If 
an enzyme is deposited into a gap of size £, the proba- 
bility of self-desorbing after a finite time is given by (see 
Eq. (Ha)) 

ndosorb(^, f)^'^-jl tanh-^ {jTJl) ■ 
Therefore, the required expression for (fr) is 

(frif)) = (21) 

Jo 

In the numerator of Eq. ([2T|) . IIdcsorb(^, /)^e~^d£ is 
the fraction of deposited enzymes that (i) land in a gap 
that has a length between £ and £ + d£ and (ii) eventu- 
ally self-desorb after finite time. In the denominator, 
/o°° ndcsorb(^, .f)^e~^d£ is the fraction of deposited en- 
zymes that self-desorb after a finite time. The result 
pTjl is confirmed by simulation data in Fig. [5^b). 

Equation (j2ip was derived by considering the deposi- 
tion of a single enzyme onto an infinite DNA with expo- 
nentially distributed gap lengths. However, as is the case 
with Eq. (|16[) . it is also approximately true for a small 
number of depositions onto a finite DNA: providing the 
number of oxoGs annihilated is small compared to the 
total number of oxoGs, the distribution of gap lengths is 



still approximately exponential. Hence, for a given de- 
position rate of enzymes per unit length onto an infinite 
DNA, Eq. (j21[l will hold approximately for times such 
that the fraction of oxoGs annihilated is small. 

Given a deposition rate of enzymes (per unit length), 
we can estimate a desorption rate (per unit length) from 
Eq. ((2T|) . If desorption were a Poisson process, then 
the desorption rate, fcoff, is found from the inverse of the 
mean unbinding time of the repair enzyme. Although 
the desorption process in our model depends on the dy- 
namics of electron charge transport (rendering it to be 
non-Poisson) , the inverse of the ensemble averaged con- 
ditional return time of an electron 1/ (f,.), is nonetheless 
a reasonable definition for the detachment rate koS- We 
expect this value of fcoff to be accurate, as long as the 
fraction of oxoGs annihilated by repair enzymes is small. 
The probabilities and times relevant to electron dynamics 
are summarized in Table [Till 



B. Colocalization of enzymes to lesions 

We now consider a permanent lesion on the DNA (one 
that does not annihilate upon absorption of an electron) . 
Such a lesion may be bound to other enzymes and cofac- 
tors so that it can act as a sink for multiple electrons, 
or it can reflect electrons. In this section, we consider le- 
sions that can cither absorb or reflect electrons, as shown 
in Fig. [TJc). We are primarily interested in the average 
number of depositions required for a repair enzyme to 
be adsorbed within a certain (small) distance from the 
lesion. 

For the sequential deposition of many enzymes onto 
a DNA populated with guanine radicals and lesions, 
the evolution of enzyme and guanine densities is not 
amenable to exact analytical solution. Therefore, our 
approach will be to track enzyme-lesion distances and en- 
zyme concentrations on the DNA by performing Monte- 
Carlo simulations. 

Each simulation consists of a series of adiabatic depo- 
sitions. A deposition is simply the spontaneous appear- 
ance of a MutY-[4Fc-4S]"^'^ enzyme at a randomly chosen 
position along the DNA. Note that a deposition is an at- 
tempted adsorption: it can result either in the enzyme 
sticking to the DNA, or desorbing from it. In our sim- 
ulations, the number of enzymes on the DNA can grow 
without bound. We do not model the bulk dynamics for 
MutY-[4Fe-4S]^"'' enzymes in solution. 

In our model, each enzyme that is deposited releases an 
electron along the DNA. However, rather than perform- 
ing time-consuming, explicit simulations of a Broadwell 
process, we exploit our analytic results to implicitly ac- 
count for the electrons. The rules for enzyme desorption 
and adsorption come from the probabilities Hdesorb and 
Hadsorb found in Eqs. pT|) . Specifically, consider the de- 
position of an enzyme, E, between two already adsorbed 
enzymes, Ei and E2 (see Fig. IDJa)). Let the distance 
from E to Ei be di, i = 1,2. Then the probability of E 
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adsorbing and knocking off Ei is ■^j^rjj- and the prob- 
ability of E self-desorbing is ^{j^Ych 



fd2 



■). In the 



case where an enzyme is deposited between a lesion and 
an adsorbed enzyme (see Fig. Eljb)), the adsorption and 
desorption probabilities have to be modified. If Ei is 
replaced by an electron-reflecting lesion, the probabil- 
ity of E permanently adsorbing without displacing E2 
is zero. The probability of E adsorbing and knocking 
oS E2 is i- 



1 



21+757 '^^'^ probability of self-dcsorption is 
1 jf ig replaced by an electron-absorbing le- 



±4- 

2^2 1+fd 
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FIG. 6: (a) Deposition of new enzyme E (solid hexagon) in 
between two adsorbed enzymes E\ and E2 (empty hexagons) . 
(b) Deposition of a new enzyme E between a lesion and an ad- 
sorbed enzyme, £2- Adsorption and desorption probabilities 
are given in Tables HVl and IVl 

sion, the probability of E permanently adsorbing without 
displacing E2 is | i+fdi ' probability of E adsorbing 
and knocking off E2 is ^ i+fd., ^'^^ probability of self- 
desorption is ^(y^j^ + These probabilities are 

summarized in Tables ITVl and fVl 



Event: 



Probability: 



E self-desorbs 



1 + fdi 



+ 



fd2 



E adsorbs. 
El desorbs 



1 



1 

2 i+fdi 



E adsorbs, 
E2 desorbs 



1 



2 l+/d2 



TABLE IV: Adsorption and desorption probabilities in Fig. 
[6l^a) when the enzyme E is deposited between enzymes Ei 
and E2. 



Event: 



Probability: 
(reflecting 
lesion) 



Probability: 
(absorbing 
lesion) 



E self-desorbs 



- + - , 

2 ^ 2 V l + /d2 



fd2 



fdi ^ fd 



1 + fdi 



l + fd2 



E adsorbs, 
E2 stays 
adsorbed 



1 



1 

2 i+fdi 



E adsorbs, 
E2 desorbs 



1 1 



2 l+/d2 



1 1 



2 l + /d2 



TABLE V: Adsorption and desorption probabilities in Fig. 
^h) when the lesion is an electron absorber and reflector. 

MC simulations were performed on a periodic domain 
of size F containing a single lesion, which is equivalent 
to a single finite domain with length F and lesions at 
X = and x = T. We start our simulations with no ad- 
sorbed BER enzyme (MutY), but with a unit density of 
guanine radicals (oxoG) whose gaps follow an exponen- 
tial distribution (see Eq. (fTSI) ). When a single enzyme 



is deposited randomly on [0,F], the positions of the two 
particles (either oxoGs, lesions or already adsorbed en- 
zymes) on either side are recorded and di and c?2 are 
calculated (see Fig. Using the probabilities in Tables 
IIVI and IVl the outcome of this deposition event is deter- 
mined: either the newly deposited enzyme adsorbs, or 
it desorbs due to its electron returning. Note that if an 
adsorption occurs, exactly one of three other events also 
has to occur: (i) a neighboring enzyme is reduced and 
desorbs (ii) a neighboring oxoG is annihilated or (iii) an 
electron is absorbed by a neighboring lesion. 

Figure [7] shows density profiles obtained from our MC 
simulations. In Fig. [TJa), the depletion of guanine radi- 
cals is greater away from lesions: a guanine radical that is 
close to a lesion can, essentially, only be annihilated from 
one side. Near x = 0, the probability of oxoGs being an- 
nihilated from the left by a rightward-moving electron is 
very small. Similarly, near x = 5, the probability that 
oxoGs are annihilated from the right by leftward-moving 
electrons is also very small. 

Figure [Tljb) shows that electron reflecting lesions even- 
tually prevent the build up of enzymes near lesions. 
The presence of an electron-reflecting lesion increases the 
local self-desorption rate. Note that the enzyme self- 
desorption probability is always greater in Fig. [6{b) than 
it is in Fig. EJa) ~ when the lesion is electron reflect- 
ing. Therefore, near a reflecting lesion, the recruitment 
of enzymes by guanine radicals has to compete with this 
increased self-desorption rate. Although the density near 
the lesion increases with time, for a fixed time, its value 
is always smaller than the bulk value. Another way to 
understand the enzyme depletion is through a particle 
conservation argument. Since the total number of gua- 
nine radicals and BER enzymes is conserved, an increase 
in oxoG density near the boundaries must correspond to 
a decrease in the enzyme density. 

Figures [TJc) and \T[^d) show density profiles near 
electron-absorbing lesions. The oxoG densities in Fig. 
[Tljc) remain essentially unchanged from those surrounded 
by electron-reflecting lesions (Fig. [Zla)). As shown in 
Figs. [Tljb) and El^d), the BER enzyme density profiles 
arc also similar for a small number of depositions, away 
from lesions. On the other hand, [7|^d) also shows that 
for larger deposition numbers, the BER enzyme density 
near electron-absorbing lesions increases markedly. 

The total number of particles on the DNA strand can 
be found by integrating the densities from a; = to a; = F. 
For example, in Fig. [7]Jb), the solid curve representing 
the enzyme density after one attempted deposition takes 
the value ~ 0.12 over most of the domain and decreases 
slightly near the lesions. Therefore the (average) number 
of enzymes that remain adsorbed after one attempted 
deposition is approximately 0.12 x 5 = 0.6. This is in 
excellent agreement with the solid curve in Fig. 3] and 
Eq. (HH) for / = 1 since (Badsorb) = eEi(l) = 0.596... 

Figure [7] only shows the densities up to 20 deposi- 
tion attempts. When the number of depositions is much 
greater than 20, all of the enzyme-seeding guanine radi- 
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FIG. 7: Implicit-electron Monte-Carlo simulations of the evolution of mean guanine radical ((a) and (c)) and BER enzyme 
((b) and (d)) density profiles after 1 (dotted), 7 (dot-dashed), 14 (dashed), and 20 (solid) enzyme depositions, (a) and (b) 
correspond to electron-reflecting lesions at a; = and a; = F = 5, and (c) and (d) are for electron-absorbing lesions. Results 
were obtained from averaging 10^ trials and using a flip rate of / = 1. 



cals arc annihilated. In the absence of any electron ab- 
sorbers on the DNA, there can be no net increase in 
enzyme number, and the enzyme density in Fig. [TUb) 
eventually saturates to unity everywhere in the domain, 
identical to the initial oxoG density. Each guanine radical 
is eventually replaced by a BER enzyme, so the long-time 
BER enzyme density mimics the initial oxoG density. 

In contrast, when the lesions are electron absorbing, 
there are always two permanent electron absorbers in the 
system. In this case, the number of enzymes can grow 
without bound, even when all the oxoGs are depleted. 

Figure [SJa) shows how enzymes converge to electron 
absorbing lesions located at a; = and .t = F = 5. At 
any given time, we label the m enzymes on the DNA 
according to their position Ei so that < i?i < £'2 < 
... < i?m < F. Both the number of enzymes on the 
DNA, m, and their positions, are functions of n, 
the number of (attempted) depositions that have oc- 
curred. We plot the quantities x\ = min(£'i,F — 
X2 = min(£'2,r - Em-x) and X3 = min(£'3,F - Ern-2) 
as functions of deposition number n in Fig. ^sl). When 
fewer than 3 enzymes are adsorbed on the DNA, we de- 
fine Xi = F, i = 1, 2, 3. From our simulations, we find the 
scaling 

~ for i 1,2,3, (22) 

in the large n limit. For a BER enzyme to successfully 
excise a lesion, we assume that it has to be within a few 



base pairs of it. We set the physical enzyme- lesion dis- 
tance Xi = xi/p ^ 5a, where a is the width of a base pair 
which we take to be 0.34 nm, and estimate n. Approxi- 
mately 1 in 40,000 guanine bases are guanine radicals^^ 
so p = (160, 000a)~^, and the number of attempted de- 
positions required for the closest sticking enzyme to be 
within 5 base pairs of the lesion is n « 6 x 10®. If each 
deposition takes at least 0.0005 seconds?^ this amounts 
to a total (minimum) search time of about 50 minutes. 
Although this is a significant reduction compared to the 
original ID sliding search time discussed in the Intro- 
duction, it is likely that MutY locates lesions even more 
quickly through a combination of the CT mechanism and 
facilitated diffusion along the DNA strand. 

The solid curve in Fig. [5]Jb) shows the total number 
of enzymes on the DNA as a function of the deposition 
number when the lesions at a; = and a:; = 5 are elec- 
tron absorbing. Upon depletion of the guanine radicals 
(shown by the dotted curve dropping to < 10~^), the en- 
zyme number increases as 

0(ni/3). The dashed curve m 
Fig. [HJb) shows the number of electrons absorbed by the 
lesion. Initially, this is less than the enzyme number since 
enzymes adsorb mainly by oxoG annihilation. However, 
as all the radicals are used up, the dashed curve asymp- 
totes to the curve for the enzyme total, indicating that 
the net increase of enzymes on the DNA is due primarily 
to lesion-induced colocalization. 
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FIG. 8: (a) Convergence of repair enzymes to an electron 
absorbing lesion. The distance between the lesion and the 
closest, second closest and third closest enzymes (denoted by 
Xj, j = 1, 2, 3 respectively) scales as 0{n~^^^) for n ^ 1 where 
n is the deposition number. Results were obtained using / = 1 
and by averaging over 5000 trials, (b) The total number of 
enzymes and guanine radicals on the DNA and the number 
of electrons absorbed by the lesion as a function of deposition 
number after averaging over 100 trials. The enzyme number 
scales as 0{n^^^) for n S> 1. The dimensionless flip rate was 
/ = 1, and the domain size was F = 5. There were initially 5 
guanine radicals present. 



Given that the enzyme-lesion distance scales as 
0{n~^/^) for electron absorbing lesions, one can directly 
show that the number of enzymes on the DNA scales as 
0{n^^^) through a simple argument. If the enzyme-lesion 
distance is 0{n~^^^) it takes 0(n^/^) attempts before an 
added enzyme lands closest to the lesion. When n is 
large, the enzymc-lcsion distance is small and the elec- 
tron released by the newly deposited enzyme will be ab- 
sorbed. For every 0{n^^^) depositions, on average, one 
permanent adsorption occurs. Hence, for every 0{n) de- 
positions, 0{n^/^) adsorptions occur. 

While the convergence of CT enzymes towards lesions 
scales as Xi ^ n~^/^, the convergence of passive enzymes 
(those that simply adsorb onto DNA without emitting 
electrons) scales . The faster convergence of 

passive enzymes^ is a consequence of linearly increasing 
the passive enzyme density on the DNA. The CT mecha- 
nism on the other hand, prevents the recruitment of large 
numbers of BER enzyme on the DNA at any given time 
with the total number scaling as 0{'n}/^) <^ 0{n) for 
large n. Hence, although BER enzymes only colocalize 
near lesions (note the maxima in the enzyme density oc- 
cur at the lesions in Fig. [7fd)), and the CT mechanism 
suppresses the wasteful build-up of enzymes in undam- 
aged parts of the DNA. 



IV. SUMMARY AND CONCLUSIONS 

We developed a mathematical model for a proposed 
charge-transport mediated mechanism of Base Excision 



Repair (BER) enzyme colocalization to DNA lesions. En- 
zymes adsorb and desorb through a charge transport 
(CT) mechanism^i^ which we model using a stochastic 
Broadwell process. Our main finding is that the CT 
mechanism concentrates repair enzymes at lesions pro- 
vided the lesions are electron absorbing. 

We first calculated enzyme sticking probabilities and 
self-desorption rates in the absence of lesions. Our results 
for an infinite, lesion free DNA, populated with guanine 
radicals, are summarized in Figs. [4] (which predicts the 
enzyme sticking probability) and (which predicts the 
electron's mean conditional return time). For the de- 
position of a single enzyme onto an infinite DNA, the 
results are exact; for a given deposition rate per unit 
length, we expect the results to hold approximately pro- 
viding the fraction of guanine radicals (oxoGs) annihi- 
lated is small. We also explored how enzymes colocalize 
to lesions using Monte-Carlo simulations. Enzymes were 
adiabatically deposited onto a circular DNA with a single 
lesion. We found that electron-absorbing lesions colocal- 
ize CT enzymes, and while electron-reflecting lesions do 
not (Fig. O. 

Simple faciliated diffusion is often unable to account 
for the fast search times observed in certain DNA-protein 
reactionsii^ Cherstvy et ali^ state that under realistic 
conditions, facilitated diffusion cannot occur and propose 
that acceleration is achieved through the collective be- 
havior of proteins. In the context of target search by en- 
zymes, the CT mechanism complements facilitated diffu- 
sion modelsiiiii^iii The CT-mediated mechanism is one 
such example of collective behavior. Instead of basing the 
enzyme search problem on the time for a single protein 
to find its target, the CT mechanism relies on a collec- 
tive build-up of enzyme density at the lesion. Hence, 
issues important in facilitated diffusion theories, such as 
the starting point of the enzyme relative to the lesion 
and the length of the DNA become irrelevant in the CT 
mechanism. 

In the case where targets (lesions) are electron absorb- 
ing, we find that the maximum enzyme density always 
occurs at the permanent lesions and furthermore that 
the CT mechanism maintains a low density of enzymes 
far from lesions to suppress oxoGs, which are another 
form (albeit less permanent) of DNA damage. In fact, 
after an initial transient where all oxoGs are annihilated, 
the density of enzymes for most of the DNA will be of 
the order of the oxoG density, which is very low (about 
1 in 160,000 base pairs). Subsequent enzyme depositions 
will colocalize only near the lesion. Our results show that 
although n ^ 10^ (attempted) depositions are required 
for the concentration to build up to a sufficient level at 
the lesions in order for them to be excised, the number of 
enzymes actually adsorbed on the DNA is much less, at 
Cl(ni/3) « 100. Although this is a significant reduction, 
it is still greater than the copy number of MutY (^ 20), 
so it is likely that the effects of ID diffusion of MutY are 
important.— 

When considering the collective behavior of enzymes. 
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one important constraint is that the number of BER 
enzymes available to participate in the search mecha- 
nism is fixed. The copy number for MutY, in particu- 
lar, is about 20^ placing a bound on the total num- 
ber of enzymes that can be successfully adsorbed on the 
DNA strand. Thus, the CT search mechanism is effective 
only if the number of oxoGs is not significantly greater 
than ~ 20. Although guanine radicals absorb electrons, 
thereby seeding the adsorption of BER enzymes, too 
many radicals can deplete the reservoir of BER enzyme 
before they significantly concentrate to the lesions. 

Although in our model, there are two modes of en- 
zyme recruitment - oxoG-mediated and lesion-mediated 
(when the lesion is electron absorbing) - it is the lat- 
ter that colocalizes enzymes to lesions. We re-emphasize 
that the initial recruitment by guanine radicals can only 
increase the enzyme density to a level that is of the order 
of the initial radical density. This density is far too low 
to ensure reliable excision of the lesion. However, upon 
subsequent depositions, enzymes rapidly colocalize and 
the accumulation is more focused. 

Although our simple model successfully predicts colo- 
calization of CT BER enzymes to electron-absorbing 
DNA lesions, it neglects many potentially important as- 
pects. For example, BER enzymes are not point particles 
but have a finite size of about 10-15 base pairs. Random 
adsorption of finite sized particles has been studied^ and 
could be used to enhance our current model. We also ne- 
glected the sliding of BER enzymes on DNA. Inclusion 
of finite size effects and enzyme sliding into our model is 
likely to decrease the search time to a lesion. The effect 
of other proteins on the DNA, besides BER enzymes, is 
also important. These proteins could physically prevent 
the adsorption of BER enzymes, absorb electrons emitted 
by BER enzymes or shield the lesion from electrons (or 
possibly all three). We currently do not know the effect 
of molecular crowding on the CT model, but this topic is 
discussed by Li et ali^ One possible approach to study- 
ing these more subtle attributes is to develop and analyze 
them within coarse-grained, mass-action type models, in 
conjunction with Monte-Carlo simulations. 
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APPENDIX A: SOLUTION OF THE ONE-SIDED 
BROADWELL PROBLEM 

Taking the Laplace transform of Eq. ^ , we obtain 

d{x - xo) 





i9Q(x, s) 
dx 



MQ(2;,s) 



(Al) 



where Q{x,s) = {Q+{x, s),Q-{x, s))'^ , Q±{x,s) 
/o°° Q±{X: t)e^''*dt and 



M 



-is + ^i + f) f 



-f 



(A2) 



The solution to Eq. (|A1[) . can be found in two separate 
regions x > xq and x < xq and matching the solutions 
with the appropriate jump conditions derived from inte- 
grating Eq. (|Aip over an infinitessimal segment centered 
about xq: 



Q+{xo,s) - Q+{xq ,s) =1, 
Q-{x^,s) - Q-{xq ,s) =0. 



(A3) 



The general solution of Eq. (|A1[) . Q{x,s;xo,£), can be 
expressed in the form 



Q(x,s) = < 




Cl 



C2 
1 

C2 



where Ai_2(s), ci,2(s) are given by 



Ai,2(s) = ±v/(s + Ai)(s + M + 2/), 
ci,2(s) 



X < Xo, 
X > Xo, 

(A4) 
(A5) 



s+A'+/-Ai.2(s) ■ 



The constants j4>, By, A<, and _B< are obtained by 
imposing the Laplace Transformed boundary conditions 
Q+{Q,s) = Q-{£,s) = 0, which come from Eq. ([9]), and 
the jump conditions (|A3|) : 



A, 



A- 



ciC2e 



-(Ai+A2)2-o rpAif+A22;o 



(Ci - C2)(cie^i^ - 026^2^) 

l^^f,^Q-{>^l+>^2)xo^QXlXO+\2i _ gAl^ + A22;oJ 

(ci - C2)(cie^i^ - C2e^2f) 

A22;ol 



(A6) 



C2e^='^[c2e"^i^° 



By = 



(ci-C2)(cie^i^-C2e^^0 ' 
(ci - C2)(cie^i^ - cae^^f) 
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APPENDIX B: LIMITING CASES OF THE 
BROADWELL MODEL 

Upon eliminating P_ from Eqs. ([l]), satisfies 



QT2 



-2{F + M) 



dP+ 
dT 



M^P+. (Bl) 



Similarly, eliminating P+ from Eqs. ([T]) gives Eq. (jBip 
but with P+ replaced with P_ . Upon neglecting electron 
decay, M = 0, and Eq. (|Bip simplifies to 



2/ 



dP+ 



d^P+ 



(B2) 



where we have used the nondimensionalization ([5]) and 
the non-dimensional flip rate / = F/{pV). When / » 1, 
we neglect the first term in Eq. (|B2[) to obtain a diffu- 
sion equation with diffusivity 1/(2/). When / <C 1, we 
neglect the second term to obtain a wave equation with 
unit wave speed. These limits correspond to a diffusive 
and ballistic electron motion respectively. 



APPENDIX C: ADIABATIC APPROXIMATION 

Since our stochastic analysis does not account for 
electron-electron interactions, we assume "adiabatic" de- 
position of BER enzymes. An adiabatic deposition of en- 
zymes occurs when each enzyme is deposited sufficiently 
slowly so that the emitted electron completes its motion 
before the deposition of the next enzyme. At any given 
time, there is at most one traveling electron on the DNA. 

Consider Figure [HI two enzymes are deposited on ci- 
ther side of a guanine radical with the left enzyme further 
away. For this example, assume that the electrons are 
always emitted toward the radical. In an adiabatic de- 
position, the deposition of the right enzyme occurs after 
the oxoG is annihilated. The final configuration consists 
of an adsorbed right enzyme and a desorbed left enzyme. 
In a non-adiabatic deposition, the right enzyme can be 
deposited before the annihilation of the oxoG. The fi- 
nal enzyme configuration depends critically on the time 
between the first and second depositions. If this time is 
long (a "late" second deposition), the oxoG is annihilated 
by the rightward electron and the final configuration is 
identical to the adiabatic case. If the inter-deposition 
time is short (an "early" second deposition) , the leftward 
electron can annihilate the oxoG first and the final con- 
figuration corresponds to an adsorbed left enzyme and a 
desorbed right enzyme. 

For a deposition to be adiabatic, the electron dynamics 
must be much faster than that of enzyme depositions: 



P 



(CI) 



where p is the density of guanine radicals and k^n is an 
intrinsic enzyme deposition rate per unit length of DNA. 



Thus, the adiabatic limit arises when p^V/kon oo and 
pF/kon OO, with / = F/[pV) fixed (to keep the over- 
all probabilities Hadsorb , Hdosorb unchanged in Eq. (fTT|) ). 
Note that / can still be small in an adiabatic deposition, 
as is the case in Fig. [9] 




Non-adiabatic 
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from left 
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deposition 
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deposition 

<} ^ 
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from right 



T 

O 



enzyme 
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-<> 



FIG. 9: Possible outcomes from an adiabatic and non- 
adiabatic deposition of a pair of repair enzymes. The left 
enzyme is always deposited first, but is further away from the 
guanine radical than the right one. The flip rate F satisfles 
F <C pV so that the electron motion is ballistic. The final 
configuration of a non-adiabatic deposition depends critically 
on the time between the first and second depositions. 



APPENDIX D: GUANINE GAP DISTRIBUTION 

Consider a lattice made up of n sites on which guanine 
radicals can randomly appear at a rate of 57 radicals per 
unit time T, per lattice site. Each lattice site can hold 
at most one guanine radical. The size of the gap be- 
tween two guanine radicals is the number of empty sites 
between them. Let N{m, T) denote the total number of 
gaps of size m (measured in lattice sites) at time T . Then 
7V(m, T) obeysS 



1 dN{m,T) 

n df 



2 iV(m',T) 

m'— m+1 



mN{m,T). (Dl) 



We will take the continuum limit of Eq. (|D1[) when the 
number of sites becomes infinite, the guanine radicals 
become points on a line and the gap length becomes a 
continuous random variable, taking any value between 
and oo. We aim to calculate the probability distribution 
function (PDF) of the gap length given a fixed average 
density of guanine radicals p. 

Let Lq be the total length of the lattice and a be the 
width of a single lattice site so that Lq = na. Further- 
more, the time taken for G guanines to appear on the 
lattice is To, where nflTa ~ G and p = G/ Lq. 

Now wc define dimensionless variables y, t and p = 



14 



p{y^ t) where 



y = pam, 
t = T/To, 

p = N/{nnT) = N/{Gt). 



(D2) 
(D3) 
(D4) 



Note that < y < oo and that for large G, Gt is approx- 
imately the total number of gaps at time i; hence p in 
Eq. (|D4p is the fraction of gaps that have size N at time 
t. 

The desired continuum limit is now obtained by taking 
n oo, ap — > so that y in Eq. (|D2[) becomes a con- 
tinuous variable ranging from to oo. and G,Lo ^ oo: 
the number of radicals that appear and the DNA length 
become infinite in such a way that p = G/Lq stays a 
constant. When these limits are taken, p{y,t) becomes 
the probability of finding a gap of length y at time t 
and J^p{y,t)dy = 1. Upon setting q{y,t) = tp{y,t), we 
obtain the intcgro-differential equation 



dt 



q{y',t)dy' - yq{y,t). 



The Laplace transform in t of Eq. (|D5P is 

/•oo 

sqiy,s)^2 q{y',s)dy' -yq{y,s), (D6) 



where q{y,s) = e 
initial condition p{y, 0) 
with respect to y gives 

dq{y,s. 



^^q{y, t)dt and we have used the 
0. Differentiating Eq. ([Pej) 



dy 



3q{y, s) 
(y + s) 



-0, 



(D7) 



which is solved by q{y, s) = A[s)/{y -\- s)^, To determine 
the integration constant A{s), we take the y — > limit of 
Eq. (jPS]) to obtain 



dq 
'dt 



q{y',t)dy' =.2t, 



(D8) 



where the last equality arises from the normalization of 
p{y, t). The Laplace transform of Eq. (|D8|) gives q{Q, s) = 
Hence, q{y,s) = 2/{y + s)^ resulting in q{y,t) = 
(D5) t^e"*** and p{y,t) = te~^*. Therefore, if Y is the non- 



dimcnsionlized gap length at t = 1, we find 



Prob(i/ < Y < y + dy) = e ^dy. 



(D9) 
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For E. Coli, the maximum deposition rate can be ob- 
tained by assuming a nucleoid radius of approximately 



h « 0.3/.im. Upon assuming a MutY diffusivity of 73 ~ 
3 X 10~^cm^/s, the Debye-Smoluchowski estimate is fcon ~ 
47rD6 f« 6 X 10^° M"^s"^ For MutY concentration of 
C f» 20 enzymes/fL, the average time between depositions 
is {kC)~^ ^ 0.0005s. 

While passive enzymes converge more quickly to enzymes 
when measured in terms of the deposition number n, CT 
repair enzymes converge more quickly when measured in 
terms of the number of adsorptions m = 0{n^^^). In fact, 
the enzyme-lesion distance for repair enzymes scales as 
0{n~'^'^) = 0(m-^) compared to 0(m-^) for passive en- 
zymes with m = n (every deposition results in an adsorp- 
tion). 



